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bigrating has previously been developed and implemented. A 
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I. INTRODUCTION 


A. BACKGROUND 

Interest in the enhanced absorption of light by a rough 
metallic surface can be traced back to 1902 when Wood first 
noted abrupt losses of approximately 90% in the intensity of 
light reflected from a metallic grating for a narrow band of 
frequencies at a given angle of incidence [Ref. l:p. 661]. 
These reflectance 'dips' were unexplained by contemporary 
theories and came to be known as Wood's anomalies. In his 
paper to the Journal of the Optical Society of America in 
1941 [Ref. 2], Fano first postulated the theory that the 
loss in the energy of the beam during reflection was due to 
‘the generation of what he termed "polarized quasi-stationary 
waves" which propagated along the surface of the metal. In 
1976 Maystre and Petit presented theoretical arguments that 
total absorption of the incident electromagnetic energy was 
possible with metallic gratings [Ref. 3]. Hutley and 
Maystre presented experimental evidence of such total 
absorption Later that same year [Ref. 4]. 

Concurrent with the more recent work in the absorption 
of electromagnetic energy by metallic gratings, there has 
been a great deal of interest and investigation of surface 
electromagnetic enhancement phenomena such as surface 


enhanced Raman scattering and enhanced second harmonic 
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generation [Ref. S:p. 1240; Ref. 62 pp. 355-367 | ee eo oes 
field enhancements on the order of 1O- have been noted in 
the literature [Ref. S5:pp. 1244-1247]. The increases in 
field enhancement occur under the same conditions as the 
increases in absorption for metallic gratings [Ref. 5:p. 
12a 

This phenomenon of energy absorption and electromagnetic 
field enhancement is due mainly to the coupling of the 
incident electromagnetic waves into surface electromagnetic 
waves in the electron plasma of the metal in the form of 
surface plasmon polaritons propagating parallel to the 
interface of the media [Ref. l:p. 683; Ref. 6:pp. 362-368). 
A surface electromagnetic wave is characterized by the 
exponential decay of its associated fields in the directions 
normal to the surface of the medium. Since the electron 
plasma on the surface of a metal constitutes a polarizable 
medium against the background of the lattice ions, an 
electromagnetic wave incident on the surface will induce a 
polarization in the medium and will in turn be modified by 
the polarization?. In such a medium, this coupled 
excitation mode is termed a surface plasmon polariton 
[Ref. 7:p. 1]. Under certain conditions of grating surface 


profile and incident frequency and polarization geometry, 


the details of this mechanism and the role played by 
the grating in causing resonance coupling are explained more 
fully in the following chapter on theory. 
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Eneweoupwng 1S resonant and a significant fraction of the 
incident energy is transferred to surface plasmon 
polaritons. This interpretation of the electromagnetic 
MimeraeeLONs ate the grating Surface has much theoretical and 
experimental support [Ref. 8:pp. 69-73]. 

Besides its obvious use in surface science in the study 
of surface electromagnetic interactions, the capability of 
the grating to couple an incident bulk light wave into a 
surface wave has given rise to a host of practical 
applications [Ref. 8:pp. 74-75]. Grating couplers can be 
used to couple a light wave into a surface wave or a guided 
wave in electro-optical devices. The dependence of the 
resonant coupling condition on the frequency of the incident 
light and on the periodicity of the grating corrugations 
also make grating couplers useful elements in filters and 
reflectors for solid state lasers and other integrated 
optics applications. These couplers can couple a laser beam 
into and out of an optical waveguide (as a guided wave) 
(Ref. l:pp. 701-702]. This capability could conceivably be 
exploited for high efficiency coupling in optical 
communications systems, some of which are presently under 
development for use as military communications systems. The 
enhancement effects available with gratings find application 
in surface catalysis in electrochemistry, in advanced 
dielectric breakdown due to surface roughness in high energy 


lasers, and in the production of radiation from a 


als: 


non-relativistic charged particle beam interacting with a 
surface grating [Ref. 8:pp. 72-73; Rete Jap. 3,) koe eee, 
Ber. 1408 

Until the middle 1970's investigations in this area 
focused almost exclusively on classical gratings, i.e. 
gratings with periodic corrugations along a single 
dimension. The most significant absorption and enhancement 
effects have been observed using gratings of highly 
conducting materials such as gold, silver and copper with 
periods on the order of the incident wavelength and 
corrugation depths on the order of one to fifty nanometers 
[Ref. 8:pp. 88-89}. These effects are highly dependent on 
the orientation of the incident electromagnetic field. For 
total absorption in classical metallic gratings, it is 
necessary for the plane of incidence to be perpendicular to 
the grating corrugations and for the incident wave to be 
linearly polarized with the magnetic field vector transverse 
to the plane of incidence. 

The dependence of the grating coupler efficiency upon 
the polarization and orientation of the incident light can 
theoretically be reduced through the Userciebaigrar mgs see 
gratings with periodic corrugations along Ssoctnh sureage 
dimensions {Ref 9:p. 2}. Within the last decade there has 
been a significant amount of theoretical work with 
bigratings [Ref..12:pp...227, 275-276, 279] to scompamisem 


to that accomplished for classical gratings, experimental 
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investigations of absorption and enhancement using 
bigratings is as yet somewhat seareee 

Although the degree of coupling is generally very 
sensitive to the angle of incidence, even in bigratings, the 
restriction on the incident polarization would effectively 
be removed. hems of bigratings, then, would allow many 
applications to be pursued in which the requirement for 
linear polarization would be impractical, such as the 
absorption and storage of sunlight as an energy source 


pPRee. 9p. 2]. 


B. PROBLEM STATEMENT 

Analysis of electromagnetic interactions at the 
interface of a grating with an air or vacuum medium is 
extremely complicated due to the geometry of the grating. 
The analysis is performed for the purpose of investigating 
surface excitations and for designing gratings. The 
approaches taken to perform the analysis nonperturbatively 
have taken two basic forms, differential and integral 
fRef. 13:pp. 15-40]. Both use the complex dielectric 
function as a frequency dependent quantity allowing 
extension of the theory to anv dielectric medium. The 
differential approaches numerically integrate Maxwell's 
equations across the grating surface boundary. The integral 


approaches use the Rayleigh method or Green's theorem. 


ey 


The one common factor in all of these approaches is 
that, though they give results which accurately conform to 
the observed experimental results, they are extremely 
lengthy and can only be accomplished through the use of 
considerable amounts of computer resources [Ref. 9:pp. 3-4]. 
The nonperturbative analysis of the bigrating is even more 
complex than for the classical grating due to the added 
dimension of the boundary conditions. For this reason 
several schemes have been proposed for perturbative analysis 
of the interface problem for a classical grating to provide 
more readily realizable results. Among these are the 
methods developed by Kroger and Kretschmann [Ref. 14], Toigo 
et al. [Ref. 15], Mills [Ref. 16]),oBeeoneandeSungaifet. aa 
and Glass, Weber and Mills [Ref. 18]. 

The perturbative technique advanced by Glass, Weber and 
Mills used an approximation to first order in the surface 
profile amplitude. The perturbation results for the 
dispersion curve were compared favorably to that obtained 
using an integral method of nonperturbative analysis 
employed by Toigo et al. Glass, Maradudin and Celli 
extended the nonperturbative analysis method for application 
to bigratings [Ref. 5]. Glass and Maradudin [Ref. 19] 
compared results from this nonperturbawiyews-eanigucweemei— 
experimental results of Inagaki et al. [Ref. 20] for total 
absorption on a bigrating with favorable conclusions. Glass 


(Ref. 21] modified and extended the perturbation technique 
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of Glass, Weber and Mills for use with a bigrating, allowing 
for arbitrary polarization, plane of incidence, and 
Simultaneous stimulation of two surface polaritons 
propagating in non-collinear directions. The results from 
this perturbation theory for reflectance and enhancment 
versus angle of incidence were compared with results from 
the nonperturbative analysis technique of Glass, Maradudin, 
and Celli for bigratings. The results compared favorably 
within the limits of the approximations upon which the 
theory was based“. 

Glass further modified his perturbation technique to 
allow for the simultaneous excitation of four surface 
polaritons and to increase its facility in treating cases of 
Horna Incidemce on a bigrating [Ref. 9]. The 
generalization also allows for determination of the complex 
dispersion relation for surface polaritons in cases where 
there may be wavevectors at the intersection of two 
Brillouin zone boundaries. The implementation of the latest 
perturbation theory advanced by Glass and comparison with 
results from the nonperturbative analysis technique of 


Glass, Maradudin and Celli is the present concern. 


The approximations and their limits are discussed in 
more detail in the following chapter on theory. 


Lo 


To ee EEO 


A. MECHANISMS 

The conduction electrons near the surface of a good 
conductor, such as silver or gold, may be considered as a 
plasma of nearly free electrons [Ref. 22:pp. 160-161]. 
Wave-like longitudinal oscillations in the charge density of 
the electron gas exist as normal modes of the classical 
system or as elementary excitations in the quantum 
mechanical system [Ref. 22:pp. 260-262; Ref. 23:pp. 
200-202]. ~Sucn quanta, CONSIStING Of 42 ecoupmed wea heme am: 
charge oscillation and electric wave, are termed bulk 
plasmons, or simply plasmons. When localized to a surface, 
as an evanescent wave, the charge density oscillation, with 
its accompanying macroscopic polarization and 
electromagnetic wave, has a component transverse to the 
direction of propagation (and hence, transverse to the 
surface). sucha SUFrace eClectromadmerle wae) coup mecmine 
collective oscillations of the conduction electrons is 
called a surface plasmon polariton. Henceforth, the use of 
the term surface polariton will be unde@rstocd) ce mean 
surface plasmon polariton. 

Surface electromagnetic waves may be stimulated by 
incident electromagnetic radiation in the form of a bulk 


light wave only in the presence of a grating. Consideration 
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Hebe ls eresenreued towthe case of a linearly polarized 
incident light wave. The periodicity of the grating surface 
profile defines a reciprocal lattice [Ref. 2l:p 2648]. As 
a consequence of the Bloch theorem, a grating reciprocal 
lattice vector may add to the component of the incident 
wavevector parallel to the surface to equal the wavevector 
of the evanescent wave at a given frequency [Ref. 22:pp. 
163-164; Ref. 24:pp. 2-26; Ref. S:p. 1241]. The resonance 
condition is acheived when the frequency of the evanescent 
wave and that of the incident light wave are equal. The 
amount of energy transferred to the surface polariton 


through this coupling can reach significant proportions. 


B. ANALYTICAL METHODS 

Since the wavelength of the surface polariton in the 
optical region is much greater than the Fermi wavelength for 
the materials of interest, the surface waves may be treated 
by classical, macroscopic, electromagnetic theory. 

The term exact theory will be used henceforth in lieu of 
the term nonperturbative analysis. The development and the 
equations used in both the exact theory and in the 
oerturbation theory are taken directly from the paper bv 
Picea ve varaducin, and Ceili (Retr. 5])and, to a greater 
extent, from the paper [Ref. 21] and the technical report 
(Ref. 9] by Glass. The exact theory is well summarized in 


Maradudin's review article [Ref. 25:pp. 423-469]. Although 


a 


the developments in both of the following subsections 
(II.B.1. and II.B.2.) should be taken as summaries of the 
analyses given in these papers, the references will not be 
explicitly cited within the subsections to avoid excessive 
repetition. 

The overwhelming majority of the literature in this area 
employs the Gaussian, or CGS, SyStemvorvuni cosas bars 
convention will be adhered to for uniformity and comparison. 

1. Exact Theory 

A schematic representation of the physical geometry 
is shown in Figure 1. The coordinate system, shown elevated 


in the figure for clarity, is lecated with ther x7 = 0 


3 
plane as the average position of the surface in the vertical 
direction. The Ky and KX, axeS are coincident with the 
orthogonal dimensions of periodicity of the bigrating 
surface. The surface profile is defined by 


=> > _ a 
x = X 


The region above the bigrating, X > Cer), is a vacuum and 
the region X. < Ck) is the dielectric characterized by 
the frequency dependent complex dielectric function 

é(w) = Ep (o) el E,(w)- 
The periodicity of the) bigqrating 1S Gesen sete eeen aeeres 


“~ a 


> ~ 

vectors a, = 4, * and a, = a, X, where ay and a, are 
the periods of the grating corrugations along each of the 
two surface dimensions. The surface profile function is 


periodic in two directions and one may therefore write 


aa 





x 


~{ 
\ 
@ 


—_— 
as 


® 


ee 
wa 





\ 


Le Ss 
AS. ‘ la Mj 
SRR 
DRSOAERES STEERS 
TESS OYPRES IFRS 
SO to 7 
ERSTE NS ER 
OOTP STF REELS 
TEES RESELL BRL 
we LBL OSL RRR a TTF. 
BRR RPE 
LESSIG BR TIPS) 
SPIRIT RR TRS 


ALPES 





N 
PRY le 









PRQuULCC Jin oeMemaulc OF Bigrating 
and Incidence Geometry. 


25 


~4 > > 
C(x) = (xX, +n a, 


+m a), 
where n and m are integers. The restriction to the specific 
case of the square bigrating, where a, = 4, = 4, will be 
imposed later. 

Light of angular frequency w is incident on the 
bigrating from the vacuum with wavevector kK at an angle 0 


from the normal. The plane of incidence is rotated from the 


axis x, by the azimuthal angle @ and contains the unit 


i 
normal a The projection of the incident wavevector Kk onto 
the X= O plane is given by 
Rk, = 2 [sine cose x ind sine x l 
ae [sin8 co X, + Sino sin Xo], (is) 
since |K| = =. 
e 


It is convenient at this point to define some terms 
for use in describing the polarization of the incident 
electromagnetic wave, in referring to specific regions in 
the area of the grating surface, and in characterizing the 
grating surface. The case of the linearly polarized 
incident wave with the electric field vector perpendicular 
or transverse to the plane of incidence (TE) will be denoted 
by the term s polarized. The case of the linearly polarized 
incident wave with the magnetic field vector transverse to 
the plane of incidence (TM) will be denoted by the term 
p polarized. The ratio of the maximum Netone of tne Suriace 
profile above the x, = 0 plane, ¢ 


3 
the grating corrugations, a, (in either direction for the 


, to the period of 
max 
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square bigrating) is defined as the corrugation strength of 
the grating. The region between the maximum and minimum 

heights of the surface profile, ¢ 
¢ = =e 


The cornerstone of the analysis is the use of the 


< X3 ¢¢ where 


nee Max 


is termed the selvedge region. 


min max ’ 


Rayleigh hypothesis. The assumption here is that 
expressions for the fields that are valid above the selvedge 
region in the vacuum and those that are valid below the 
selvedge region in the dielectric medium may be extended 
into the surface itself. This assumption places a 
limitation on the validity of the analysis and this 
limitation applies to the perturbative analysis as well. If 
the corrugation strength becomes large enough, field 
components backscattered from the sides of the corrugations 
down into the corrugation troughs become significant and the 
analysis breaks down. As noted by Glass, Maradudin and 
Celli, however, the hypothesis has been shown to be valid 
outside its normal limits for periodic surface profiles 
which are analytic. 

In the Rayleigh method, an exact expression is 
written for the electric field above the selvedge region 
which satisfies Maxwell's equations in the vacuum and 
See mies ene Silochy condition for the doubly periodic 
geometry of the bigrating. This expression is written as a 
Rayleigh expansion in terms of the incident and scattered 


fields. Through the use of the vectorial equivalent of the 


Zo 


Kirchoff Integral for diffraction, coupled with the 
extinction theorem, the expression for the field in the 
dielectric medium below the selvedge region is eliminated by 
rewriting it in terms of that above the selvedge region. 
This halves the number of unknowns in the problem and is 
termed the reduced Rayleigh method. The boundary conditions 
associated with Maxwell's equations are then applied at the 
interface using the expression for the fields above the 
selvedge region, which is an implicit application of the 
Rayleigh hypothesis. The results are two sets of linear 
algebraic equations, of infinite extent, with the Rayleigh 
coefficients as the unknowns. The Rayleigh coefficients 
determine the magnitudes of the scattered field components 
for both the diffracted waves and the evanescent waves. 

To solve these equations numerically, the matrices 
must of course be truncated to some degree. The procedure 
in the calculations is to truncate the equations to a given 
finite dimension and solve them numerically. The equations 
are then truncated to some larger finite dimension and the 
numerical solution is repeated. Iterations of this 
procedure are carried out until convergence of results is 
either confirmed, in which case the validity of the Rayleigh 
nypothesis for the particular case being investigated is 
taken to be established, or until divergence becomes 
apparent. In the latter instance some other analysis 


technique must be attempted in order to investigate the case 
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of interest. For cases where convergence iS apparent, a 
Suitable dimension for the matrix equations is chosen 
dependent upon the degree of accuracy required. 

The quantitative analysis is initiated by writing 
the electric field in the vacuum above the selvedge region 
in the form 
H(w,X) = B*(0,R,) exp[ik, ° xX 


+ bs E’(w,Ka) exp[ika ° K, + tao(w,Ka) x3]. (2) 


The summation in the second term is over all the translation 
weceors)) G, of the reciprocal to the lattice defined by the 
geometry of the square bigrating. The reciprocal lattice 
vector, then, can be obtained with the expression 


ma 2 a a 
—— (mM) xX, 3° MX 


a TOmeeiiowe= O77 +], -t2,.«. «—. C3) 


2) 
In equation (2), the wavevectors for the incident wave, k, 
and for each of the scattered waves, Ra, are used in the 


forms given by 


k = R, - ag(o-k,), (4a) 
and | 
Ra = Ka + ao(w,Ka), (4b) 
where 
Ko =. eee (4c) 


G I 
BAewQuanticy Aq (wo, Ke) is then given by the expression 
fo? 7c? - Kal’, fire KA < ical (5a) 
Aq (w, Ka) 7 ; 
i [x - annem HO) KA > BOVE, GSipy 


Zar 


and the quantity Ao (w,k,) may be obtained by evaluating 
equation, (5) for, Uneseanewune ine G=0 and Ka = i, : 


ao(W,Ky) = ao (w/Ka) — ~ cos@ . (6) 


The quantities Rt and E> in equation (2) denote the vector 
amplitudes of the incident and scattered fields, 
respectively. The vector amplitude of the incident field is 


given by 


k 
it 7? + . a : 
E (w,K, ) = i, + X, — By + ey x a Bie (7) 


where BY and By determine the magnitudes of the p and s 
components Of me incidemt ficeldjeseseoeeively ene ey cece 
amplitude of the scattered field is given by 


be a K 
EB (0 /Ka) = Rg = ss A, (Ka) 


“3 Go (0, Ka) 

+ BE ¢ Kal A, (o,Ka), (89 
where the A, and AY determine the magnitudes of the p and s 
components of the scattered field, respectively. These are 
the Rayleigh coefficients which are the unknowns in the 
problem. Each term in the summation in equation (2), i. e. 
each scattered field, represents either a diffracted beam, 
when Ka < = , or an evanescent wave, when Ke > = : 

Glass, Maradudin and Celli are quitesdecailed jin 

their method of quantitative elimination of the expression 


for the field within the dielectric medium and derivation of 


the final set of linear equations for the Rayleigh 


Zo 


coefficients. The vectorial Kirchoff integral provides an 
effective boundary condition for the field in the vacuum. 
AS an application of the Rayleigh hypothesis, the expression 
for the field above the selvedge region, equation (2), is 
then used in the boundary condition equation. 

At this point Fourier expansions in terms of the 
Surface profile are introduced with the intention of 
rewriting the vector integral as a set of linear equations. 
One of these expansions is critical to the development of 


the perturbation theory and is therefore stated here: 


exp[-ial (X,) ] = Dn ¢(a|G) expiG e se (9) 


In tnis expansion, # is defined by the integral expression 


9 (a|G) = a {| dx, dx., exp|-i{2 e ae “ at (2 )}], CO) 
C a 
S 


where an is the area of a unit cell of the bigrating 
surface. The quantity a inside the dielectric must be 
distinguished from the quantity ag in the vacuum, due to the 
complex dielectric function, e(w), of the medium, and is 


given by 
a(w,K2) = je(w) — sz - K : (les) 
c a2 e 


CoS se Natecucineandeesl limause. che, Femrier 
expansions to write the Kirchoff integral equation as a 
doubly infinite set of simultaneous, linear, inhomogeneous 
equations for A, (Ka) and A, (w, Ka). Using the 


abbreviations of A, (G). FOr A, (o Ka) and A, (G) for A, (o,/Ka), 
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the linear equations are given by 


Az, C - CG! 
a thee 2 2) aaa. A, (G’) + baa, a, (87) 
fe 
— seg 1 3} laa, B, + bay B, |. (12a) 


ae 


$ [a A Gu 3°] | 
a Se aaa, cag. A, (G’) - aga, a, (8) | 
__ {eg 8p 


- Ba lea By - ao B, |. (1 Zon) 


with the following definitions of terms: 


eeu = a(w,Ke) - Bo (0 ,Ka,) , (i3a)} 
Ba = a(w,Ke) + Ao (w,k, ) P (1315s) 
aae. = K we Ee x Ka]. (lag) 
baa, = Ka & Ka, ; (13a) 

Ke Ke, 
Ca@- = Daa, + a(@,Ka) ao(0,/Ka,) ’ ae 

and 

a —— ll 13f 
Cann i a(o,Ka) ao (0-k,) a 


2. Perturbation Theory 


The essence of the» perturbation Eseery is co seduce 
all expressions to first order results in terms of the 
surface profile ¢. This process is begun by expanding the 


integral #¥ to first order in ¢( yielding the result 
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g(a|@) x SG gis vaGcc), (eae 


where 


eC \e= -= | dx, dx, C (X,) exp|-i Cue z, | (15) 
e a 
Cc 


The use of the Kronecker delta function is introduced into 
the notation in equation (14). With the assumption that the 
dimensionless parameter af is small, the expression for # 
given by equation (14) may be used in the equations (12) for 
the Rayleigh coefficients. With the further assumption that 
0.) = 0, the resulting first order approximations for A 


Il 
and A, are given by 


[ee ° i: a 
A, (G) = et ba(t-8a,e [Fe cagA, (2 ) - aae,A,(G + 


aaa | 5a 9 - iba C(C) 
aa ee | lea = Co B, |. eee) 
and 
A, (G) = iaae a(t-8a,a-) [FSB r{agg.ay (8 - baa-A, (2’)}] 


6a 5 - tba C(G) 
- aaa ee lagy By + Dao B,|- (86,15)) 
The assumption is then made that the Rayleigh 
coefficients of the specularly diffracted beam, A, (0) and 
A, (0), are of a magnitude such that the amplitude of the 
specular beam is much greater than any other diffracted 


beam. The second assumption made here is that there is a 


possibility of up to four evanescent waves in a condition of 


ou 


resonance corresponding to four resonantly stimulated 
surface polaritons. The amplitudes of a resonant evanescent 
wave will be much greater than that of any of the 
nonresonant terms. The wavevector for such an evanescent 


wave is given by 


Ka = R, if G, = Repo) for r= 1,2,3,4, (17) 


where Rep (@) is the wavevector of the resonantly excited 
surface polariton of frequency w, and where w is equal to 
the frequency of the incident light. 

Since treatment is limited to only small deviations 
from the flat surface, it is beneficial to discuss the 
limiting case of oe) = 0 to gain some insights to the 
physical situation. In the flat surface limit, which is 
zeroth order in the surface profile, the dispersion relation 
for the surface polariton is given by 

AO Ney + €(W) to (@/K.) = Q. (18) 


The dielectric medium is assumed to be isotropic so that to 


zeroth order in (, the dispersion relation depends only on 


the magnitude of the surface polariton wavevector, Ksp° For 
a given frequency, then, Ksp will describe a circle in 
k-space on the K, = 0 plane. Figure 2 is a schematic of 


this situation with the solid circle representingea carcic 


of radius Ksp° The wavevector, Ka , of a resonantly 
i 


3This limitation has been imposed with the use of the 
Rayleigh hypothesis. See pages 25 and 26 in the exact 
theory section for a dicussion of this limitation. 


5) 7 


excited evanescent wave will have its tip at this circle and 
will satisfy the dispersion relation given by 


a(w,Ke ) + E(w) Ao (W,Ke ) = QO. Loy 
i iE 


The present formulation allows for the possibility of four 
different reciprocal lattice vectors (Gq through G, bral 

Pigure 2) coupling to K, with the resulting evanescent waves 
satisfying equation (19). The particular situation depicted 
in Figure 2 shows a case where only one of the resulting 


evanescent waves, Ka , satisfies the resonance condition. 
4 


The limitation of this schematic representation is that the 
Meeaprocal lattice vectors for the bigrating exist only if 
Eee surrtace is not flat, in’ which case the constant 
frequency circle for Ksp is distorted ra some extent by the 
bigrating corrugations. 

Since surface polaritons on a flat surface are p 
polarized, only the A, coefficients (not the A, 
coefficients) of the resonant evanescent waves need be 
considered as large with respect to the nonresonant terms. 
Thus the assumption is that the Rayleigh coefficients A, (0), 
Bo). BCD) 7 A, (G,), A, (G3), and A (G,), may be larger 
than all others and should therefore be treated differently. 
Meese ecoeadement muse allow for the possibility that there may 
be no resonantly excited evanescent waves and the 
corresponding Rayleigh coefficients would then be of the 
same magnitude for the specially treated terms as for all 


the others. 


BS 
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Figure 2. Schematic of Wavevector Coupling. 


Having reviewed the basic premise of the 
perturbation theory, a descriptive summary of the rest of 
the development advanced by Glass will now be employed. A 
statement of the final results will be given at the end of 


this summary. 
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A further abbreviation is introduced into the 
notation at this point and mention is made of it here since 
it is used in the statement of the final results. In this 
abbreviation, G, is replaced everywhere with its subscript 


alone, so that aaa 1s now written as a ; A, (G3) as 
Lag 


eZ 
A, (3), and so forth. 

Equation (16) is used to write an explicit 
expression for each of the six important terms. Each 
important term is thus written in terms of the other five 
and also of a sum over the nonresonant terms. Each given 
nonresonant coefficient can in turn be expressed, by 
equation (16), in terms of the six important coefficients 
and all the other nonresonant coefficients, The latter are 
dropped: each given nonresonant coefficient is expressed 
only in terms of the six important coefficients. These 
equations for the nonresonant coefficients are then 
substituted into the expressions for the resonant and 
specular beam coefficients. 

Within the development of the perturbation theory, 
an equation arises for the frequency of the surface 
pemartton to zeroth order in (¢. This expression, which is 
the solution of equation (18), will be useful in later 


Giscussions and is stated here as 


2 L+e Dee, 
oa | S KE : (265 
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The final result is a set of six linear algebraic 
equations for the resonant and specular beam coefficients in 


and B of the 


II ee 
These equations are given in matrix form as 


terms of the amplitude coefficients, B 


incident beam. 


Deo. 
Cae oe “ae A a oT ey 
aD 
Le ay ae A He ND aay 
a2 
Me De) eS a ne ue eqs 
ee, 
A “42 AS iat) Ms Ny ieee 
| Ry Ro R, 5 S (1+T) | JA, (0) | 
Smee me 
Se 2) aa 
ee | ee oe 
——s (U,-M,)B, - N,B, (21) 
(WHV-P)B, - OB, 
a 
00 
(X-S)B, = ee 
It should be noted that the quantity 9, defined by 
- w ra 
0 = —gre/a = xX 22) 


is merely the incident wave frequency expressed as a 


dimensionless quantity. 


defined by 


Also note that the quantity Q is 


(23) 


where O, is just the surface polariton frequency for the 
flat surface, O- from equation (20), renormalized by the 
nonresonant terms, which were kept in the theory. An 
explicit definition for oO. is given in equation A.1 of 
Appendix A. The other terms appearing here are defined 
explicitly in Appendix A taken from the technical report by 
Glass. 

The procedure of calculation, then, is to solve this 
set of six equations numerically for the resonant and 
specular beam coefficients. These results are then used in 
the equations for the nonresonant coefficients to obtain 


their values. 
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IIIT. NUMERICAL IMPLEMENTATION 


The numerical implementations here are, for purposes of 
comparison, exactly analogous to those performed by Glass 
(Ref. 21] in instituting the perturbation theory with 
Simultaneous excitation of two resonant surface polaritons. 
In his quantitative evaluations employing the exact theory, 
Glass parallels those conducted by Glass, Maradudin and 
Celli [Ref. 5]. The exact method here duplicates that 
employed by Glass but for a somewhat different geometry. 

The implementation of the four-polariton perturbation theory 
will parallel the method used by Glass in the implementation 


of his original) perturbation tmeonr 7. 


A. QUANTITIES AND GEOMETRY INVESTIGATED 
In the numerical investigation of the theory, normalized 


values are chosen for the incident amplitude coefficients, 


By and B,- In the cases investigated, either total p 
polarization, By, = 1 and B, = 0, or tOcal s polarizariony, 
By = 0 and By = 1, was chosen for the incident beam. The 


resulting Rayleigh coefficients, then, will be some fraction 
of unity indicating the magnitudes of the p and s components 
of the evanescent waves and the diffracted beams relative to 


the magnitude of the incident beam. 


Se 


ine surtacesprofile used for the investigations 
duplicates the one employed by Glass [Ref. 21}. It 
describes a square sinusoidal bigrating with a profile 
symmetric along xy and Xo Wichiepemlodicity a. “The bigrating 
Surface is depicted schematically in Figure 1 in Chapter II. 


The surface profile function is given by 


Zt 
«| + cos| 5 x] 
tr an 
“= A, {cos | 7 ‘7 + cos | = x |} 


~ Aid {cos |= eo cos | 27 x,]} (24) 


The reciprocal lattice vectors for this square bigrating are 








=> ss 27 
C(x) ) = A, {cos| = 








given by 


G(l,m) = - f. aa oT bal L,m = OQ, BURP E 2 poss ® (25) 


Thus one may refer to G(l,m) as the (l,m) reciprocal lattice 
vector. Similarly, a particular scattered beam, 
corresponding to the (l,m) lattice vector, may be referred 
to as the (l,m) diffracted beam or the (l,m) evanescent 
wave. With this reference scheme, the wavevector of a 


scattered beam projected onto the surface is written as 


= aad 27rl ~* 21m ~*~ 
K(l,m) = K, + G(L,m) = K, + —— x, + SS x,. (26) 


Once numerical results for the Rayleigh coefficients are 
Seeained, equation (8) is used to find the amplitude of the 
scattered field. The amplitude of the incident field is 
obtained with equation (7) using the known amplitude 
coefficients B, and B,. The reflectance of the (l,m) 


Il At 
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diffracted beam [Ref. 21:p. 2653] is then given by 


2 
IBSto,R(L,m) ]| aofw,K(l,m) ] 
ee (27a) 


Z 
#024] ao [o-k, | 


with the total reflectance given by the expression 





R= d, d RCL on ae (27b) 
where the primes on the indices indicate summations are 
carried out only for l’ and m’ such that Raa < ae 
thereby excluding evanescent wave terms. Using the 
amplitudes of the scattered and incident fields, the 
amplitude for the total field may be obtained by evaluating 
equation (2). The electric field enhancement is defined 


[ReEs 21:p. 26 5seas 


JB(o,%) |7 
Ve (28) 
1 > 
Bee) | 
3 *max 
where {Siem indicates a point just above the selvedge region. 


B. EXACT THEORY IMPLEMENTATION 

The method of exact theory implementation described here 
is a summary of the procedure outlined by Glass [Ref. 21;pp. 
2653-2654]. This reference will not be cited further within 
this section. 

In determining the Rayleigh coefficient values using the 


exact analysis, the integral # given in equation (10) must 
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be evaluated for all possible reciprocal lattice vectors. 
There is, in general, no analytic expression for # and 
numerical calculations thus require performance of a fast 
Fourier transform (FFT) for each possible G. The surface 
profile function term with the height coefficient Aya: the 
cross term, would require the use of a two dimensional 
Hounrer tEoamnstorm in direct application. Expansion of the 
exponential of the cross term may be used to circumvent the 
use of a two dimensional Fourier transform. With this 


method the expression for the integral # for the cross term 


case may be rewritten as 


mn 
09 —ten 
$(a|@) = a pepe g'™) alo.) ¢™) (alG,), (29a) 


where 
G=G, x, + G, oa (29b) 
and 
+a/2 
$(alG) = - | cos” [2% exp|-ic,x| 


-a/2 
+ exp|-ta{ hy cos | 22%] + h, cos [42%] }] dx. (29) 
This integral can then be evaluated using a one dimensional 
SPT and a numerical result may be obtained for %(a|G) 
truncating the summation in equation (29a) to as many terms 
as is required to acheive the desired degree of accuracy. 


Using this numerical result for g9(a|G), equation (9) is 


then solved as a set of linear equations in matrix form in 
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order to determine the values for the Rayleigh coefficients. 
To do this, as previously indicated, the resulting doubly 
infinite matrix equation must be truncated to some degree. 
The dimensions of the matrix determine the reciprocal 
lattice vectors that are retained in the numerical 
calculation. For truncation to a matrix of dimensions 

2N* x Nn“, the lattice vectors G(l ,m) retained in the 


calculation are those where l and m satisfy the condition 


Na 


a Se tm ee 5 


5 < (30 
Thus the diffracted beams and evanescent waves corresponding 
to reciprocal lattice vectors with either index outside 
these limits is ignored in determination of numerical 
results. For the limitation of both l and m to the range 
from -2 to +2, the resulting N is 5 and the matrix to be 
solved is then of 50 x 50 dimension; for the range -3 to +3, 
N = 7 and the matrix is 98 x 98; for -4 to +4, N = 9 with a 
162 x 162 matrix to be solved. Should accuracy require 
extension of the limits for l and m to -5 to +5, N is then 
ll and the matrix dimensions are 242 x 242. 

The elements of the matrix are complex, as are the 
required FFT's. Using the full (double) precision available 
with a 32 bit machine* is desirable to obtain the necessary 


accuracy, particularly in the larger matrix calculations. 


‘the machine employed here was an IBM 3033 system and a 
32 point mesh was used in the FFT's for the exact 
calculations. 
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With such requirements in computational resources for 
numerical determinations based on the exact analysis, the 
need for reliable perturbation schemes is clearly evident. 
The Fortran code employed in calculations here is the 
same as that used by Glass in his exact calculations for 


comparison with the two-polariton perturbation theory. 


C. PERTURBATION THEORY IMPLEMENTATION 

Since the integral 9(a|@) has been eliminated from the 
calculations in the perturbative analysis, no FFT's need be 
performed. The quantity Lic defined by equation (15), is 
all that is required and, for the surface profile defined by 


equation (24), iS written [Ref. 21:p.2654] as 


Y 7 
Bae) = 5 51 [hy 2 \le|| saws 2 5 in| ,2| 


1 
as. 5a o(My aie 5111021 


+ 


alee nil 2 eee 


The set of linear equations defined by the 6 x 6 matrix 
equation (21) are then solved to determine the values for 
the Rayleigh coefficients. 

The numerical calculations here still involve complex 
quantities and Full (double) precision was also employed. 
The calculations, however, by excluding FFT's and reducing 


the linear equation solution to a 6x 6 problem, require 
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computational resources that are orders of magnitude less 
than the exact calculation requirements. 

The Fortran code employed for the perturbation 
calculations is a modification to that used by Glass for the 
two polariton perturbation code. Necessary changes were 
Carried out and the results tested against exact results for 


identical cases in order to correct program errors. 


D. SELECTION OF PARAMETERS 

To conform with the experimental work of Inagaki et al. 
[Ref. 20] and with the theoretical work of Glass and 
Maradudin [Ref. 19] and that of Glass [Ref. 21], the 
wavelength of the incident light was chosen as 633.0 nm. In 
order to test the possibility of four-fold resonance 
advanced by Glass [Ref. 9:p. 13], it was desirable to choose 
the geometry so as to have resonant absorption at normal 
incidence. In this manner, resonance effects would be 
dramatic, i. e., approximately total absorption of light at 
normal incidence, and any qualitative trends in resonance 
conditions would be readily apparent. 

For normal incidence, the projection of the incident 
wavevector is zero and the resonance condition, from 


equation (17), is then 
K =0+Cs= Repl): (32) 


Choosing the (1,0) lattice point for convenience, the 
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resulting reciprocal lattice vector from equation (25) is 
given by 


an_ey (33) 


Calor ~~ X, 





Combining these, the resulting resonance condition is 
Z 
Z Z, Z 21 
Kop = Kp = 67 (1,0) = [2] : : (34) 


Using a zeroth order approximation in ¢, equation (20) may 


be employed with the result from equation (34) yielding 


2 
oe mE. MAE c* [2] . (35) 
ig é a 


Assuming that the frequency of the incident light equals the 
frequency of the surface polariton at Ksp = 27/a, and 
eeeuming that this surface polariton frequency is 
approximately equal to the flat surface polariton frequency, 
then w can be replaced by w in equation (35). Then one may 


write 
Zz 
Zon | oe é 
a? = [225] |---|. vale) 
Substituting the wavelength of the incident light in vacuun, 


Ao, into this equation, the final result is obtained for the 


approximate periodicity for resonance at normal incidence as 


ate [—+_]*. (37) 


AS in the work of Inagaki et al., Glass and Maradudin, and 
Glass, the material chosen for the bigrating was silver. 


Interpolating from the data measured by Johnson and Christy 


[Ref. 26] for Ag, the value of the dielectric function at 
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the incident wavelength of 633.0 nm is € = - 18.3 + i 0.479. 
Taking the real part of this value the resulting value fora 
from equation (37) is 615.47 nm. 

Although the projection of the incident wavevector onto 
the surface is zero for the case of normal incidence, some 
finite non-zero value must be used for IK, | in the numerical 
calculations in order to define the unit vectors in equation 
(17). In executing a computational run, the.value for Ky is 
not defined as an explicit input parameter but is derived 
from the specification of the angle of incidence, 6. For 


2 degrees for 


this reason 9 was specified as 10. 
computational runs in cases of normal incidence. 

Taking the values of a = 614.47 nm and Ag = 633.00 nm, 
numerical calculations of exact theory were performed for p 
polarized light at normal incidence and zero azimuth with 
the surface profile height coefficients h, and hii as zero. 
This case was repeated for several choices of hy in order to 
establish an optimum value for maximum absorptance at normal 
incidence. This resulted in the choice for the optimum 
coupling value of A, as 7.4 nm. The results of these 


calculations are not formally presented here as their sole 


purvose was to establish this optimum coupling value. 
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IV. VALIDATION AT OFF-NORMAL INCIDENCE 


Within the literature there has been considerable use of 
two types of scans in the numerical search for the resonance 
SeecietOn OL a Darticular case {Ref. 27]. In the first 
method the reflectance is scanned for a resonance dip by 
Variation of the angle of incidence, 6. Referring to Figure 
Hin subsection B.1 of Chapter II, it is apparent that this 
is an implicit variation of the value of ka As the value 
of Ky is changed, the resonant coupling condition is 
approacned, met, and then passed and the reflectance is 
scanned through the resonance dip. A second method is to 
vary the incident frequency. Referring to Figure 2 in 
subsection B.2 of Chapter II, this has the effect of 
changing the radius of the constant frequency circle at the 
tip of the resonant surface polariton. Here, K, is held 
fixed and the resonant coupling condition changes, thereby 
scanning the reflectance as a function of frequency. For 
cases of off-normal incidence, only scans versus the angle 


of incidence have been executed. 


Along witn the scans for reflectance, the concomitant 


enhancement peaks were also obtained. These results are not 


"see equation (1) in subsection II.B.1. 
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central to the present discussions and are therefore 


relegated to Appendix B as Figures 28 through 46. 


A. RESULTS AT OFF-NORMAL INCIDENCE 

As discussed in section D of Chapter III, the values of 
A = 633.00 nm and a = 615.47 nm were chosen for the 
incident wavelength and bigrating Sarton and then the 
optimum coupllingevaluc on 7-4 nieron A, was established at 
normal incidence. A value of 5 degrees was arbitrarily 
chosen for the azimuthal angle, #, and scans of reflectance 
versus the angle of incidence, 9, were carried out for 
several values of A, - Tne primary purpose in these 
calculations was to establish the limits of validity for the 
perturbation theory in cases of off-normal incidence. The 
values of hy used were multiplicative factors of the 
baseline value: one-third, 2.5 nm; one-half, 3.7 nm; unity, 
7.4 mm; and twice, 14.8 nm. For all of these calculations, 
the values of h, and Ali were held at zero. From equation 
(24), the resulting values for Ce may be calculated and 
the corresponding corrugation strengths® Abeu0r O08 rag 0 127 
0.024, and 0.048, respectively. The incident light was 
specified as s pelemizged for ali@eases ay eee ny = 722 3an 
case repeated using p polarization. The results of the 


scans for reflectance for each of these cases are 


coe equations (34) and (36) in section III.D. 
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graphically depicted in Figures 3 through 7 on the following 
five pages. Within these figures, as with all of the 
figures to be presented, the results for exact theory 
calculations .are shown as a dashed curve and the results for 
perturbation calculations are shown as a solid curve. 

As can be seen from a comparison of Bre eGaeeon the 
perturbation theory is the most reliable at the weakest 
corrugation. Minimum reflectance for this case is 
approximately 0.8. Fair agreement is held between exact and 
perturbation results for increases in hy Eons. 7 AM vane 
7.4 nm as Shown in Figures 4 , 5 and 7. In the latter case, 
jac retlectance minimum falls to approximately 0.13 for s 
polarization. With the increase in hy to 14.8 nm, the 
pertubation theory breaks down and is unable to predict the 
results, as can be seen in Figure 6. 

Some definitions for use in quantitative comparison of 
exact theory and perturbation theory results need to be 
Seo etatedeak thls point. ~The departures which are to be 
quantitatively characterized are the difference in the 
magnitude of the reflectance dip and the difference in the 06 
location of the dip. The differences should be normalized 
With an appropriate quantity so that they may be considered 
as a type of percentage error. 

For the reflectance minimum, the difference in the 
maximum absorptances for exact and perturbation results will 


be used. This quantity will then be normalized with the 
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maximum absorptance of the exact theory results. To avoid 
baeroauctiom of amother quantity, this error quantity will 
be defined in terms of reflectance minima. The reflectance 
minimum for the exact theory results is denoted as &, and 
Or the perturbation results as %g. With this notation, the 


difference as a percentage, AR, is defined as 


| (1-%o) - (1-%S) | IZ - S| 


ao age  < +00% =. 7-7 — 


| _ oe x 100% . (38) 
Characterization of the difference in the 6 location of 
the reflectance minima will be accomplished through the use 
of the differences in the magnitudes of the corresponding 
Rep wavevectors. From equations (17) and (26), it can be 
readily seen that this is identical to the differences in 
the magnitudes of the corresponding K, wavevectors. This 


difference, AK. will be normalized with the calculated 


De 
value of the surface polariton wavevector magnitude. This 
value of Ksp may be taken directly from equation (34) as 
2r/a. It can be readily determined from equation (1) for Ky 


that 


4 = W ; ns : 4 
AK SD ~ Kio - Kio =e (sin6, Simon) ; (395) 


where, as before, the subscript naughts indicate values at 
minimum vetlectance and the prime indicates the perturvation 
Eewoemvevalue. Dividing by the normalization factor, K sp’ 
and using the relations in equations (36) and (37) to 


rewrite the multiplicative constant, the normalized 


5 


wavevector difference as a percentage, which will be denoted 


as Ak, can be defined as 
Ak = a |sinOd,g - sindg| x 100% , (40) 
O 


. where Ag is the wavelength of the incident light in the 
vacuum. 

Similar definitions may be obtained for the differences 
in peak enhancement magnitudes and the peak enhancement 
locations. Although presentation of these curves is 
confined to Appendix B to promote continuity, the 
quantitative comparisons will be included in tabulated 
results within the text of the chapters for completeness. 
For the percentage difference in enhancement peaks, the 
definition employed is 


heey - ect 


Now a 2 ae 100% . (41) 
Q 


The definition of the angular error, expressed as a 
percentage in wavevector form, is exactly analogous to 
equation (40) and will be denoted by Ak . for the enhancement 
values. 

The values for the percentage differences for the five 
cases at the azimuthal angle of ® = 5° are given in Table 1, 
on the following page. Note from these results that the 
reflectance from exact calculations for the AL value of 
14.8 nm is less than 0.1 with an enhancement of over 400. 
The perturbation theory, however, predicts a reflectance dip 


down to 0.800 with an enhancement of only 99.1 for this 


56 


case. The position errors are an order of magnitude greater 
than the others noted in Table 1 but are still less than 


five percent. 


TABLE 1 
PERCENTAGE DIFFERENCES OF EXACT AND PERTURBATION 


RESULTS FOR @ = 5° AZIMUTH 


h,(nm) Pol. Ko on AS Ak Aé- Ak 
DS 5 0.808 65.0 3.5% 0.10% 2.9% Gee 
ge 7 5 0.634 Po aos 0.23% 5.5% ames 
7.4 5 0.153 342. 9.2% 0.77% 2.9% 0.84% 
14.8 5 0.095 457. | Tans 3.1% 77 .% 2.0% 
TA Pp 0.981 3.35 6.4% 0.23% 10.% 0.17% 


This breakdown supports the explanation offered by Glass 
[Ref. 21:p. 2654] for the limitation of his perturbation 
theory using two resonant polaritons. Briefly, the argument 
is as follows. The basis of the perturbation theory is the 
expansion of the quantity exp[-ia{]. At resonance, the 


approximation may be made that 


2r |e | ¢ 
leCK(1om) 1 fnax * ee al Sa 
K a 


When 0. Gees. exceeds one, the expansion series does not 
converge rapidly and truncation to a few terms is no longer 
BEaPg@odecppsOximation or its value. Thus the perturbation 


theory breaks down. 


Sif 


For the hy = 14.8 nm case, evaluation of equation (42) 
yields a value of 1.3 for CU esc aus At the 7.4 nm height, the 
value is at ae ~ 0.65. Thus the perturbation theory using 
four polariton resonance appears to follow the breakdown 
pattern of the two polariton theory in cases of off-normal 
incidence. 

A comparison of the results for AL = 7.4 nm with s and p 
polarization of incident light illustrates. an important 
point. Note in Table 1 that the reflectance minimum for the 
s polarized case is 0.15, while the minimum for the p 
polarized case is only down to 0.98. Figure 8 is a 
schematic of the wavevector coupling with an azimuth of 5°. 
In this schematic., unlike Figure 2, the two dimensional 
Brillouin zone boundary is shown, represented by the solid 
Square surrounding the constant frequency circle. Also, 
this circle is pictured with the gaps which are present at 
the intersections with the zone boundary. The directions of 
the incident fields, for the cases ofsssancenepolaniaaeuen 
are shown projected onto the kK, = 0 plane. The wavevector 
of the resonant evanescent wave is represented by the dashed 
vector with the open arrowhead (not labeled). 

The surface polaricon fas Sern a eeweteme ee ree 
component normal to the surface and longitudinal field 
component parallel to the surface and in the direction of 
polariton propagation. An incident wave of s polarization 


has its electric field vector perpendicular to the plane of 
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Figure 8. Schematic of Wavevector Coupling 


incidence and, therefore, has no field component normal to 
the surface. An incident wave of p polarization, however, 


does have an electric field component normal to the surface. 


at 5 Azimuth: 


1 


The magnitude of this normal component is proportional to 


the angle of incidence. 


a” 


Se er 


Referring to Figure 8, it can be seen that with s 
polarization, the incident field is almost collinear with 
the direction of propagation of the surface polariton and 
thus with its longitudinal field component. Whereas with p 
polarization, the incident field projection is nearly 
perpendicular to the surface polariton wavevector. Since 
the angle of incidence in these cases is very shallow, 
approximately 10 to 11°, the incident field component normal 
to the surface in the p polarized case is a very small 
fraction of the total. \Thus~ehe coupisinGg Aisemarras, 
dependent upon the collinearity of the incident field 
component parallel to the surface and that of the 
longitudinal component of the surface polariton field. As 
can be seen from the schematic, this collinearity is almost 
nonexistent with an azimuth of @ = 5°. Hence, the large 
differences in the reflectance minimums for the s and p 
polarizations of incident fields. 

At an angle of incidence of 9 = 90°, the projection of 
the incident wavevector, Ki would reach the light line, 

w = Ck. This extremum cannot, of course, be reached in 
practicality but may only be approached. The azimuthal 
angle can be made Large enough so that the resonant 
condition cannot be met at any angle of incidence. In these 
cases, Ky cannot be made large enough to position any of the 
reciprocal lattice vectors such that the coupling produces a 


resonant evanescent wave. 
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In the geometry under investigation, this limiting case 
occurs at approximately 29° of azimuth. In order to 
validate the perturbation theory more fully, reflectance was 
scanned versus the angle of incidence at an angle of azimuth 
of 25° for the single case of h, = 7.4 nm. These 
calculations were performed using both s and p polarization 
for the incident light. The results from these calculations 
are presented graphically in Figures 9 and 10 on the 


following two pages. The numerical comparisons of the exact 


TABLE 2 
PERCENTAGE DIPFERENCES OF EXACT AND PERTURBATION 


RESULTS FOR © = 25° AZIMUTH 


h, (nm) Pol. Ky ec Ak Ak Aé Ak . 
7.4 s 0.739 2577 1.3% 0.45% ie 26 0.45% 
Wage! Pp O27 76 41.5 9.1% 0.03% 19 .% 0.05% 


and perturbation results are given in Table 2. Although the 
percentage errors for the reflectance minimum and 
enhancement peak are somewhat large for the case of p 
polarized incidence, the perturbation predictions have 
angular errors that are less than 0.1%. 

Figure 11 is a schematic of the wavevector coupling for 
the 25° azimuth geometry. The magnitude of the reflectance 
minima are very nearly the same for the s and p incident 


polarization cases. Unlike the 5° azimuth cases, the angle 
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Reflectance Curves for S Polarized Light 
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Figure 11. Schematic of Wavevector Coupling 
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of incidence in these cases is rather large, approximately 
60.5°, and the p polarized incident wave has a large field 
component normal to the surface. This, together with the 
increased collinearity of its parallel component with the 


longitudinal surface polariton field COMmpOnent, Aaeeounk=s tom 


64 


the smaller reflectance minimum with p polarized incidence 
at an azimuth of 25° vice 5°. The s polarized incident 
wave, totally dependent for resonant coupling on 
collinearity with the longitudinal surface polariton field 
component, has a significantly larger reflectance minimum at 
25° azimuth than at 5° azimuth. Hence, the rather small 
differences in reflectance minimums for the s and p incident 


polarizations in the ® = 25° case. 


B. CONVERGENCE OF EXACT THEORY CALCULATIONS 

As discussed in section B of Chapter III, the numerical 
implementation of the exact theory requires truncation of 
the infinite series of equations (9) and (29a) in 
calculating the values of the Rayleigh coefficients and 
$(a|G). If the height coefficient in the surface profile 
function for the cross term, Aiy: is sufficiently small, 
then retention of a single term of the series in equation 
(29a) will suffice. In all of the cases presented within 
this chapter, the cross term is absent, nia = 0, and 
truncation of this series was not a factor in convergence of 
ene calculations. 

Convergence of the calculations with the exact theory 
here, then, is solely dependent upon the validity of the 
Rayleigh Hypothesis for the cases treated. The validity of 
the Rayleigh Hypothesis is, in turn, dependent upon the 


corrugation strength. Convergence checks were conducted for 
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the maximum corrugation strengths used in each case of 
polarizauLlonvancwagluuen combination to ensure thes aibidmay 
of the exact calculations. 

As discussed in section B of Chapter III, truncation of 
the doubly infinite set of linear equations from equation 
(9) results in elimination of reciprocal lattice vectors 
outside a chosen range of indices. Recalling the notational 
definitions from this discussion, the truncation, with the 
linear equations in matrix form, results in a matrix of 
dimensions 2N7 x 2N7 and the reciprocal lattice vectors 
retained have indices given by equation (30). 


All results presented are for N = 7. Convergence checks 


were conducted by using N as 5, 7, 9, and 11 in successive 


TABLE 3 
SUMMARY OF CONVERGENCE CHECK RESULTS FOR 
REFLECTANCES AT OFF-NORMAL INCIDENCE 


h, (am) Pol. @ Ak AR AR 


5,7 7,9 9,11 

14.8 s 5° 4.9x10° 3.9 x 10° 1,0 © log 
7.4 é 5° oe x Om 2.1 x 10° 0x10” 
Tae s Zor 5.2 x 10° 2.9 x 10/7 1.0 x 107° 
Tee D 25° 52 oes 3) Oo 2.0 x 1077 


calculations of the resonant case, i. e., at the reflectance 
minimum. Table 3 is a summary of the differences in the 


minimum reflectances. The subscripts on the labels indicate 
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the successive values of N used in the calculations. Since 
the reflectance values are already normalized with respect 
to unity, the accuracies acheived may be read directly from 
the table. From the values in the table it can be seen that 
truncation to a matrix of dimensions 98 x 98 with N = 7 
gives accuracies to at least one part in two thousand for 


the reflectance minimum. 


TABLE 4 
SUMMARY OF CONVERGENCE CHECK RESULTS FOR 


ENHANCEMENTS AT OFF-NORMAL INCIDENCE 


h,(nm) Pol. ¢ Age AZ, aby 44 tony 
14.8 S 5° 13.723 1.846 0.308 439.4 
7.4 p 5° 2.4x10°% 2.1x107? 1.8x1074 35350 
7.4 s 25° 0.339 2.8x10* 2.4x1077 a5 65 
Fea p 25° 0.175 2.3x10- 2.0xlor” 39.18 


Table 4 is a summary of the differences in enhancement 
peaks in the convergence checks. Since these are not 
normalized differences, the peak enhancement value for the 
N = 11 calculation is included in the table for each case. 
Using this value for a normalization factor, the accuracies 
acheived for the enhancements with N = 7 are at least to one 


part in fifteen hundred. 
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V. VALIDATION NEAR NORMAL INCIDENCE 


In comparing the results of perturbation theory 
calculations around normal incidence with those of the exact 
theory, several cases of incidence geometry were 
investigated. Scans of reflectance versus angle of 
incidence were performed for various surface profiles in 
attempting to obtain agreement between perturbation and 
exact results. Although the perturbation theory was unable 
to predict both reflectance and enhancement versus angle of 
incidence with an acceptable degree of accuracy for any 
Single case, some utility in predicting bigrating 
efficiencies was established for limited cases of surface 
profile and incidence geometry. Additionally, a scan versus 
incident frequency was performed for a single bigrating 


profile. 


A. RESULTS NEAR NORMAL INCIDENCE VERSUS ANGLE OF INCIDENCE 
As before, the values of A = 633.00 nm and a = 615.47 nm 

were held fixed and investigations around normal incidence 

were commenced using the optimum coupling value of 

A, = 7.4 nm with A, and Aly both at zero. Since the 

perturbation theory gave reliable results at off-normal 


incidence for this geometry, it was expected that some 


degree of reliability could be expected near normal 
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incidence as a benchmark for validation. This was not found 
to be the case. For this reason, changes in calculation 
parameters for the perturbation theory which might improve 


agreement with exact calculations were sought. 


In all previous perturbation calculations, the four 
Rayleigh coefficients which were taken to be the most 
Significant evanescent wave amplitude coefficients were 
Byt,0), A, (0,1), A, (-1,0), and A,(0,-1). This set of 
evanescent wave Rayleigh coefficients will henceforth be 
referred to as the canonical coefficients. These 
coefficients, along with the specular beam coefficients 
A, (0,0) and A, (0,0), constitute the six important terms 
referred to in the development of the perturbation theory in 
subsection B.2 of Chapter II. Inspections of printouts for 
the Rayleigh coefficients from exact theory calculations 
revealed that several evanescent amplitudes other than the 
canonical coefficients were of significant magnitude. In 
order to obtain better prediction behaviour from the 
pertubation theory, the choice of which four evanescent wave 
coefficients are included in the six important terms was 
modified. The choices for the four terms included was 
dependent on the particular case of incidence geometry. The 
Specific choices made for these terms, which will henceforth 
be termed the corrected choice of coefficients, will be 


discussed later. 
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Figures 12 through 15 on the following four pages are 
graphical representations of the reflectance versus angle of 
incidence for the perturbation and exact theory calculations 
for s and p polarization at azimuths of ® = 0° and ® = 90°. 
Table 5 below is a quantitative comparison of these results 
using the same parameters for comparison as developed in the 


preceeding chapter. In this table, the different choices 


TABLE 5 
PERCENTAGE DIFFERENCES OF EXACT AND PERTURBATION 


RESULTS FOR hy = 7.4 nm NEAR NORMAL INCIDENCE 


6 Pol. Pert. Re ac AR Ak Aé AK, 
0° s canon. 0.003 40.4 -53% 3.4% 34.% Yaege 
Os Scone 0.003 40.4 1.3% 3.1% 11.% sea 
0° p canon. 050027 3 11.% 0.17% 68.% “Ole 
O- ip “corm 0.050 273. 47.8 990502%  =SOmte 0.092 
90° s canon. 0.003 "40.4 .63% 2.5% 6.4% 2.6% 
90° “s “eer 0-003 40.4 Io cuore ll. ~~ aa 
90° p canon. 03050) 2732 122%. 0.15%. 60.20 enor One 
90°." bp Peorr. 0.050 273. 47.% 0.02% 50.%  ‘OmO9% 


ror the evanescent wave Rayleigh coefficients are indicated 
in the column headed Pert., for perturbation, as either 
canonical or corrected. From this table and from Figures 12 
through 15, it can be seen that the perturbation theory is 


unable to reliably predict results for this bigrating 
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geometry. In particular it should be noted in the figures 
that the magnitude of the reflectance dip at zero degrees 
(8 = 0.00001°) incidence is quite different in the 
perturbation and the exact theory curves. Although the 
Prgers in the angular location of the dip expressed in 
wavevector form are all less than 4%, reliance on the 
perturbation theory as a vredictor of bigrating performance 
at normal incidence could lead to misconceptions of the 
pattern of bigrating efficiency around normal incidence. 

In discussing the corrected choices of the evanescent 
wave Rayleigh coefficients, it is helpful to once again 
Consmeaer a schematic of the wavevector coupling. Figure 16 
is such a schematic for the case of near normal incidence 
with s polarization. The wavevectors of the two resonantly 
stimulated evanescent waves are again depicted as dashed 
vectors without labels. As the angle of incidence is varied 
near @ = 0°, the wavevector ki + G@ can be thought to move 
across the gap induced by the Brioullin zone boundary. The 
resonant condition is met when the reciprocal lattice 
vectors touch the constant frequency circle. At an azimuth 
of @ = 0° with s polarization, there is no incident electric 
c27-1a Component wm the (1,0) and (-1,0) directions. Thus, 
the surface waves at G, and ee are not directly excited. 
Moreover, there is no first order coupling of the excited 
waves at G, and G, to the degenerate waves at G, and G, 


since the coupling coefficients €(+1,+1) (proportional to 
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4,,) are all zero for this surface profile. On the other 


hand, the excited waves at Cs and er, do "Gouptle? row rirs = 


order, through the reciprocal lattice vectors G. a 
1,2,3,4), to other waves not on the constant frequency 


Circle. Thus a logical choice for the corrected set of 
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evanescent wave Rayleigh coefficients for this case is 

A, (0,1), A, (0,-1), A, (0,2), and A, (0,-2). With s 
polarization at an azimuth of @ = 90°, the incident field 
vector lies wholly along the (1,0), (-1,0) axis and the 
choice for the corrected coefficients in this case is 
A,(1,0), A 


(-1,0), A, (2,0), and A,,(-2,0). 


II il 

It is obvious that the two cases of @ = 0° and ® = 930° 
azimuth are identical for either polarization. With s 
polarization, for example, the incidence geometry pictured 
in Figure 16 is merely rotated 90° for the latter case. 
Indeed, a close comparison of the exact results for the two 
azimuth cases for either polarization showed that the 
reflectance curves are identical. There are, however, 
slight differences in perturbation theory results for each 
of the two azimuth cases for both s and p polarization. It 
is for this reason that the two azimuth cases were examined. 
Validity of the exact calculations has thus been reaffirmed 
and a further test of the perturbation theory has been 
acheived. 

Figure 17 on the following page is a schematic of the 
wavevector coupling with p polarized incidence. At this 
point it is necessary to note that the resonance condition 
is not sharp. As can be seen in all the reflectance curves, 
there is a degree of broadening around the theoretical 
matching condition. The breaks in the constant frequency 


Circle are narrow enough in this geometry that the resonant 


a, 





Figure 17. Schematic of Wavevector Coupling for 
P Polarization Near Normal Incidence. 


Z : : i ; ee 
reflectance dips overlap across the Seats aoe seers 


overlap that gives rise to the reflectance dips in the p 


ue discussed earlier, these breaks in the constant 
frequency cirlce are caused by the bigrating lattice and 
occur at the intersections with the Brioullin zone boundary. 
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polarized geometry at normal incidence. The exact curves 
for the s polarized cases graphically indicate this overlap 
Phanecterisic. Referring to the schematic for the s 
polarized case in Figure 16, as the angle of incidence is 
increased from zero, the Brioullin zone boundary is scanned 
from the center of the break in the constant frequency 
circle (9 = 0° on Figures 12 and 14), where the tails of the 
polaritons strongly overlap, to the edge of the break 
(9 = 1.5° on Figures 12 and 14), where the resonance 
condition is fully met. Thus, there is very low reflectance 
at @ = 0° but the minimum reflectance occurs at 9 = 1.5%. 
Returning to the ee aitauale im Figure 17 fer p polarized 
incidence, as sie angle of incidence is increased from zero, 
the scan passes right through the center of the break in the 
constant frequency circle and thus crosses the center of the 
polariton tail overlap. Therefore there is a minimum 
reflectance condition at 9 = 0° in Figures 13 and 15. Note 
that the reflectance at 6 = 0° for s polarization (Figures 
12 and 14) is identical to that at 6 = 0° for p polarization 
HatlGaves 13 and 15): both correspond to exactly the same 
excitation at the center of the break in the constant 
paegueacyechrcle. 

At ® = 0° azimuth with p polarization, the incident 
eleceric fieldwvector liceswalong the (1,0), (-1,0) axis of 
the lattice. Consequently, there is no first order coupling 


of the incident wave into a surface polariton with the 


co 


Comm and the G(0,-1) reciprocal lattice vectors. A 
corrected choice for the set of evanescent wave Rayleigh 
coefficients in this case, then, is identical to that for 
the @ = 90°, s polarized case, i. e., Bone 0) ee O el 
A, (2,0) and A, (-2,0). Similarly, the corrected choice for 
the set of coefficients for ® = 90° azimuth with p 
polarization is identical to that for the 6 = 0°, s 
polarized case. 

The preceding discussion has served to illustrate much 
of the physics involved in the coupling mechanisms for cases 
near normal incidence. However, as can be seen from the 
reflectance curve comparisons and from Gable 5, ehereuas 
very Little improvement in the perturbation theory 
predictions with the corrected choice of coefficients over 
those with the canonical coefficients. 

There are three basic assumptions in the foundation of 
the perturbation theory. The first is the validity of the 
Rayleigh hypothesis in the particular cases being 
investigated. This assumptionwis also imherent in thewexace 
theory and is taken to be valid on the strength of the rapid 
convergence of the exact calculations. The second 
assumption is that a first order truncation of he vexpansicn 
for exp[-iat] yields a valid approximation for the 
expansion. As discussed in the previous chapter, this 
assumption is taken to be valid for values of af Gal de 


max 
evidenced by the results for calculations at off-normal 
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ineidence. The third assumption of the first order 
perturbation theory is that the most important evanescent 
waves in the Bloch sum for the surface polariton are the 
four which are, on the flat surface, degenerate and whose 
wavevectors are connected by reciprocal lattice vectors; 
i. e., those evanescent terms whose wavevectors are at the 
intersection of two Brillouin zone boundaries. In the 
bigrating profiles considered up to this point, the values 
ices h, and hii have been held at zero. The reason for the 
poor performance of the perturbation theory near normal 
incidence with such a surface profile, even when al ae is 
less than one, is due to the fact that these terms are not 
present in the surface profile. For these cases, the third 
premise of the perturbation theory is not valid, since the 
coupling coefficients actually vanish among the four 
evanescent waves considered as important. Thus, there are 
other terms in the Bloch sum that are of the same order of 
magnitude as these four. In these cases, the gaps in the 
dispersion curves, at the Brillouin zone boundaries at kK = 
et, G., Gai, and G,, and the breaks in the constant frequency 
circle are due to second order coupling of the four zone 


7 7 3 er = 
coundary waves’. Thus, the use of a first order 


Sthe breaks in the constant frequency circle, depicted 
schematically in the wavevector coupling figures, are 
responsible for the shift of the reflectance minimum from 
9 = 0° to @ = 1.5° in Figures 12 and 14. 
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approximation for treating the resonant excitation of modes 
near the gap at cm Cu, cm and G, in such cases is 
inappropriate. An increase in the hoa cross term, then, 
would allow first order coupling between these four terms. 
In this manner, the magnitudes of the four evanescent waves 
considered as important in the perturbation theory, and 
which are retained in all portions of the perturbation 
calculations, could be increased over those of all the other 
terms in the Bloch wave sum. Furthermore, the gaps at G,, 
G,, Clay and G, could become first order effects. Thus, by 
increasing the cross term, the basic premises of the 
perturbation theory might once again be valid. 

Cases of surface profiles with the hy cross term 
increased from zero to 2.5 nm were investigated for 
comparisons between exact and perturbation results. The h, 
term was decreased to 6.0 nm to avoid making at ax greater 
than one, and calculations for s and p polarized incidence 
with ® = 0° azimuth were executed. The reflectance curves 
for these calculations are presented in Figures 18 through 
23 on the following six pages. Table 6 on the seventh page 
following is a listing of the quantitative comparisons of 
the perturbation and exact results for these cases. 

A definite trend of improvement for the perturbation 
theory predictions can be noted in these results. 
Particularly in a comparison of the results for the surface 


profiles of h, = 4.5 nm with and without the cross term 
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Figure 18. Reflectance Curves for hy = 6.0 nm 
and Ald = 2.5 nm with S Polarization. 
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TABLE 6 
PERCENTAGE DIFFERENCES OF EXACT AND PERTURBATION 


RESULTS FOR CROSS TERM RUNS AT © = 0° AZIMUTH 


he Leeman <0 ee Ag Ak Ag Ak, 

6.0nm 2.5nm s 0.057 37.6 .40% 1.6% 146.% 2.1% 
6.0nm 2.S5nm p 0.113 276. 11.% .00%  104.% 31% 
4.5nm 2.5nm s 0.235 29.7 2.7% .92% 28.%  .56% 
4.5mm 2.5nm op 0.265 192. 4.8%  .00% 42.% 17% 
4.5nm 0.0mm s 0.245 28.9 4.3% 1.8% 31.% 1.9% 
4.5mm 0.0nm p 0.254 179. 7.6% .00% 56.%  .03% 


present. Though further investigations were not carried out 
in this area, one may surmise that the perturbation theory 
would yield more accurate predictions for surface profiles 
in which the cross term was an increasingly important 


PraAcCker. 


B. RESULTS AT NORMAL INCIDENCE VERSUS FREQUENCY 
The surface profile for which the reflectance versus 
angle of incidence near normal incidence were predicted most 


Eoneaohyeey ele perclurodtion —Eneory, t. e., Nn, = 4.5 nm, 


—_ 


1,5 =a@mampweand 14 = 2. Was USCC tO test tne 
perturbation theory predictions for reflectance and 
enhancement versus incident frequency. Equivalently, this 
is a scan of the reflectance and enhancement versus incident 


photon energy, hw. The reflectance curves for these 
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calculations are presented in Figures 24 and 25 on the 
following two pages. Table 7 is a listing of the 
quantitative comparisons of the perturbation and exact 


theory results. The quantites used here for the reflectance 


TABLE 7 
PERCENTAGE DIFFERENCES OF EXACT AND PERTURBATION 


RESULTS VERSUS INCIDENT PHOTON ENERGY 


hy bij Pol. &5 Eo A AE Aé AE, 
4am Wee orl Ss 0.234 296 2.0% -02% 31.% -04% 
4.5nm 2.5nm Pp 0.234 29.8 1.4% -01% 14.% -O01% 


and enhancement percentage differences are the same as those 
used in comparisons of the results versus angle of 
incidence. The differences for the incident photon energies 
at which minimum reflectance and maximum enhancement occur 
are expressed as simple percentage errors. The percentage 
error in the energies for the reflectance minima are denoted 
by AE and for the enhancement peaks by AE, - From this table 
and from Figures 24 and 25, it can be seen that the 
perturbation theory results versus incident photon energy 
are appreciably better than tnose versus yenGl eee me oe mee 
for the normal incidence cases. 

The results versus incident photon energy are somewhat 


Surprising in that the dispersion curve gap does not appear. 
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Since there are gaps in the dispersion curve at the 
PeL_oullin wenesboundariesmem the we@ratinewlattice, i. e. at 
k = nw/a, n = £1,£2,..., one would expect to see a pair of 
dips in the reflectance at such a boundary rather than the 
single dip in evidence in Figures 24 and 25. In the cases 
investigated here, n = 2 and the gap concerned lies at the 
sceena Brioullin zone boundary. The reason for this 
apparent single dip in reflectance at the second Brioullin 
zone boundary is that the energy gap in the cases 
investigated is narrower than the width of the states on 
either side of the gap; i. e., the overlap of the 


neighboring states has completely bridged the gap. 


C. CONVERGENCE OF EXACT THEORY CALCULATIONS 

For the cases treated here, the maximum value of the 
cross term was Ai = 2.5 nm and this value is sufficiently 
small that the use of a single term in equation (29a) was 
effective in obtaining accurate results. 

The truncation of the infinite series of equation (9) to 
the value of N = 7 was all that was required for excellent 
degrees of accuracy”. Convergence checks were performed 
Pomme ne Maclinuim COLZUGAaLLON Strengths for each polarizations 


Although the profile with A, = 6.0 nm and Aly =2.5 nm has a 


2See section IV.B on convergence in off-normal 
calculations for a detailed discussion of the truncations 
implied here and their effect on accuracy of calculations. 


ele. 


slightly smaller corrugation strength than the hy = 7.4 nm, 
Ayy = 0 nm profile, convergence checks were conducted for 
the s and p polarizations of the profile with the non-zero 
cross term to ensure validity of the calculations. 

Table 8, below, is a summary of the differences in the 
reflectance minima for the convergence checks. As with the 


off-normal convergence checks, the differences are 


subscripted with the two N values used in the calculations. 


TABLE 8 
SUMMARY OF CONVERGENCE CHECK RESULTS FOR 
REFLECTANCES NEAR NORMAL INCIDENCE 


h, (nm) h,, (nm) Pol. A& Ak A® 


5,7 7,9 9,11 

7a 0.0 s S05°x 10° 2.1 x 10° 25x lone 
7.4 0.0 p 3.0 x 10°92 2.1 x 1072 1.3 x 107? 

6.0 2.5 s 9.3k 10> @ieeme 10° Same cm. 

eno 5 p  2.4x10% 1.6x10° 1.0x 107 


From the values listed in the table, it can be seen that the 
use of the 98 x 98 matrix in exact calculations resulted in 
accuracies to at least one part in forty-five hundred and in 
one Case CO" Getter enan Oneweane in on. 

Table 9, on the following page, is a summary of the 
differences in enhancements for the convergence checks. As 
with the off-normal convergence checks, the accuracy figure 
is obtained by normalizing the differences using the 
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SUMMARY OF CONVERGENCE CHECK RESULTS FOR 
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ENHANCEMENTS NEAR NORMAL INCIDENCE 
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eee 


Poor 
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9.5x10_ 
Se oxo 
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9.4x10- 
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oT 
eres 0 
 Beoxto 
2 8.0x1077 
eT TT O 


im Calculation. 


cy 


40.241 


oo Ore 
Seo 


Sere Oc 


The 


accuracies acheived here are at least to one part in three 


hundred fifty. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


The perturbation theory has been validated to yield 
reliable results at off-normal incidence for corrugation 
strengths of 0.024 and less. Its performance in predicting 
bigrating absorptance efficiencies and enhancement 
capabilities is limited by the quantity aE ae which must be 
less than one. Thus, Glass's modification to his 
perturbation theory to include simultaneous excitation of 
four resonant surface polaritons, rather than only two, has 
not altered its performance in cases of off-normal 
incidence. This may be verified by comparison of the 
results contained in this work with those presented by Glass 
[Ref. 21] for the two resonant surface polariton 
perturbation theory. 

Use of the perturbation theory at normal incidence will 
yield results for reflectance and enhancement versus angle 
of incidence to at least within an order of magnitude of the 
exact theory values. The degree of accuracy acheived is 
dependent upon the amount of cross coupling present for the 
particular surface profile being investigated. From the 
results presented here, the perturbation theory appears to 
perform best when the surface profile induces appreciable 
degrees of coupling through cross terms. It is recommended 


that further investigations be conducted in this area by 


AC 


performing comparison calculations using the perturbative 
and nonperturbative analysis techniques with surface profile 
functions of more evenly distributed values for the height 
coefficients. Based on the results obtained in this work, 
it is believed that a marked improvement in the performance 
of the perturbation theory would be acheived for such 
Surface profiles. There has been considerable work using 
SaerOoeEn type diffraction grating profiles [Ref. 28]. A 
bigrating with this type of surface profile may well be 
treated using the perturbation theory with much greater 
accuracy near normal incidence than has been found with the 
surface profiles dealt with here. 

In frequency scans at normal incidence the perturbation 
theory may be employed with considerably more confidence. 
Reflectance magnitudes may be relied upon to be within 
approximately 5% of the exact theory values and the 
predictions for optimum coupling frequency can be taken as 
accurate to within approximately 2%. The enhancement 
magnitudes in such scans, however, may still only be taken 
as accurate to within an order of magnitude. 

A final observation is made here that the reflectance 
vaiue for a given surface profile at exactly normal 
incidence on a bigrating is independent of the incident 
polarization. Figures 26 and 27, on the following two 
pages, illustrate this fact graphically for the surface 


profile of A, = 4.5 nm and Ai = 2.5 nm. Note in Figure 26 
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Polarizations Versus Incident Photon Energy. 


go 


that the magnitude of the reflectance curves for s and p 
polarization are identical at @ = 0°. The curves for the 
reflectance of the two incident polarizations in Figure 27 
are identical for all incident photon energies. The data 
from other calculations also show chee the magnitudes of the 
reflectance minima for exactly normal incidence are 
numerically identical for s and p polarizations to the full 


accuracies of the exact calculations. 


APPENDIX A 


DEFINITIONS OF PERTURBATION MATRIX ELEMENTS 


Equation (21), which is stated as the final result of 
the perturbation theory (subsection II.B.2), contains 
several elements whose definitions are not given in the 
text. The equations which define these terms are quoted 
here (taken from Glass [Ref. 9]). The prime on the 
summations indicate the specified sum is taken over all j 
except the specular (j = 0) and resonant wave (j = 1,2,3, 


and 4) terms. 
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APPENDIX B 


ENHANCEMENT CURVES FOR CASES INVESTIGATED 


This appendix is a collection of the curves for 
enhancement versus angle of incidence and versus incident 
photon energy. The figures are are presented here in the 
Same order in which the corresponding reflectance curves 


appear within the text of Chapters IV and V. 
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Figure 28. Enhancement Curves for 
5° Azimuth with h, at 2.5 nm. 
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Figure 30. Enhancement Curves for S Polarized Light 
at 5° Azimuth and h, at 7.4 nm. 
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Figure 31. Enhancement Curves for 
S° Azimuth with h, at 14.8 nm. 
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Figure 32. Enhancement Curves for P Polarized Light 


at 5° Azimuth and ah at 7.4 nm. 
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Figure 33. Enhancement Curves for S Polarized Light 
at 25° Azimuth and h, at 7.4 nm. 
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Pigmre 34: Enhancement Curves for P Polarized Light 
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Figure 35. Enhancement Curves for Ay =e7.4 nm Near 
Normal Incidence with 0° Azimuth and S Polarization. 


112 


Enhancement 





ak O 
: b = O 
wee i P-Polarized 
. : nN = 633nm 
Ht h, = 7.4nm 
§ 
22 | hy = O-Oim 
4 
: 1 A. = Oenm 
i 11 
200.0 + ! ' 
' 4 
| 
175.0 ! ' 
g 
' Perturbation 
‘ Canonical 
‘ aoe 
§ 
150.0 \ Perturbation 
7 Corrected 
, 
§ 
125-0 © 
§ 
§ 
§ 
i 
100.0 4! 
a, 
t 
§ 
t 
t 
75.0 
_§ 
a 
50.0 i uN 
25.0 es 
a 
O.Q beveteie tie beakslaiess oie a el terenere ise pe Mise 


GUmOnmy YO2 § OS QUO4N vols Gols 07 #108 
Angle of Incidence (@, degrees) 


Hieire@so. Enhancement Curves for h, ——aenm Near 
Normal Incidence with 0° Azimuth and P Polarization. 


five 


Enhancement 


50.0 pr SRE anRaEE RAR onan a pe 








O 
6 = 90 
1 S-Polarized 
Ao aly < A = 633nm 
h, = /.4nm 
aan iy = 0.0nm 
h, = 0.0nm 
\ 
35.0 
| ‘ \ 
30.0 ; 
, " 
25.0 | : 
| \ 
| ‘ 
20.0 | \ 
8 \ 
15.0 | : 
S 
\ 
10.0 ‘\ 
Perturbation - 
Canonical oe 
Perturbation és 
5.0 _ \ 
L&orrected = 
0.0 ~ = te re tn tnt 
O 1 2 s 4. 5 6 7 
Angle of Incidence (©, degrees) 
Figure 37. Enhancement Curves for hy = 7.4 nm Near 


Normal Incidence with 90° Azimuth and S Polarization. 


114 


Enhancement 


0.0 


a 
es 


ma 6 = 90° 
= P-Polarized 
: A = 633nm 
1 h, = 7.4nm 
: h, = 0.0nm 
h = 0.0nm 


Pertu rbation 
Canonical 


Corrected 





4 
J 
1 
J 
J 
I 
I 
1 
I 
f 
1 
1 
1 
1 
1 
i 
J 
J 
I 
! 
} 
t 
; Perturbation 
I 
1 
{ 
1 
J 
1 
1 
J 
1 
I 
I 
1 
1 
J 
J 
1 
| 
I 
I 
I 
I 





oe 
os 7 SSS 
ee ee en a a ee 
0.0 0.1 0.2 0.3 0.4 OES 0.6 0.7 0.8 


Angle of Incidence (©, degrees) 


Figure 38. Enhancement Curves for A, = 7.4 nm Near 
Normal Incidence with 90° Azimuth and P Polarization. 
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